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Ionospheric  Scintillationi/TEC  and  In-Situ  Density 
Measurements  at  an  Auroral  Location 
in  the  European  Sector 


1.  INTRODUCTION 


3 

A  wide  variety  of  C  I  systems  suffer  degradation  in  performance  due  to  phase 

and  intensity  scintillations  imposed  by  the  ionospheric  irregularities  of  electron 

density.  Earlier  publications  (Aarons)1  have  documented  the  morphological  features 

of  major  scintillation  regions.  Long-term  studies  have  predominantly  utilized 

2  3 

intensity  scintillation  measurements  (Basu  and  Basu;  Basu  and  Aarons;  Aarons 

4 

et  al  ).  The  DNA  Wideband  satellite  made  possible  phase  scintillation  observations 

5 

at  the  equator  (Livingston)  and  the  auroral  oval  in  both  the  Alaskan  sector 


(Received  for  publication  12  August  1987) 
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0  7 

(Rino  and  Matthews)  and  the  North  Atlantic  sector  (Basu  et  al)  .  Since  the  Wide¬ 
band  satellite  was  sun-synchronous,  observations  were  available  only  near  local 
midnight  and  in  the  pre-noon  hours.  To  extend  the  diurnal  coverage,  phase 

g 

scintillation  measurements  utilizing  Fleetsat  (Basu  et  al  )  and  a  quasi-stationary 

9 

polar  beacon  satellite  (Basu  et  al  )  have  been  presented. 

There  is  a  great  deal  of  interest  in  understanding  the  development  of  such 
irregularities  at  high  latitudes  where  the  ionosphere  is  often  strongly  coupled  with 
the  magnetosphere.  In  such  an  environment  the  distant  magnetosphere  serves  to 
activate  different  sources  of  free  energy,  for  example,  electron  precipitation, 
field-aligned  currents,  electric  fields,  and  so  on,  that  control  the  formation  of 
ionospheric  irregularities.  No  longer  is  it  possible  to  pursue  a  study  of  the 
irregularity  development  in  the  local  ionospheric  environment  without  considering 
the  coupling  between  the  ionosphere  and  the  magnetosphere.  This  approach  is 
particularly  useful  in  extrapolating  our  knowledge  of  the  natural  ionospheric  ir¬ 
regularity  structures  to  problems  related  to  the  structuring  of  artificially  injected 
plasma  clouds  in  the  high  latitude  ionosphere. 

In  an  effort  to  study  plasma  structuring  in  the  above  context,  the  Defense  Nuclear 
Agency  (DNA)  launched,  on  27  June  1983,  the  HiLat  satellite  in  a  circular  830-km 
orbit  at  82°  inclination.  The  satellite  transmits  coherent  signals  at  137,  390,  413, 
and  436  MHz  and  the  phase  reference  signal  at  1239  MHz  to  measure  complex  signal 
scintillation  and  total  electron  content  (TEC).  It  also  carries  a  variety  of  in-situ 
probes  providing  measurements  of  ion  density,  ion  drift,  energetic  electron  pre¬ 
cipitation,  field-aligned  currents,  and  emissions  at  two  visible  wavelengths.  All 
instruments  except  the  Langmuir  probe,  the  vacuum -ultraviolet  imager  and  a  part 
of  the  magnetometer  continue  to  operate  reliably.  Table  1  [reproduced  from 
Fremouw  and  Wittwer9 10]  lists  the  instruments  and  the  experimenters. 


6.  Rino,  C.  L.  ,  and  Matthews,  S.  J.  (1980)  On  the  morphology  of  auroral  zone 

radio  wave  scintillation,  J.  Geophys.  Res.,  85 -.4 139. 
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Table  1.  HiLat  Experiments 
(Reproduced  from  Fremouw  and  Wittwer1^) 


Instrument 

Observations 

Operating  Characteristics 

Experimenters 

Coherent 

Scintillation 

Phase -coherent  radio -wave 

Carlson  (Air  Force 

beacon 

Total  electron  content 

signals  at  138,  390,  413, 

436  and  1239  MHz; 
circular  polarization. 

Geophysics  Laboratory) 
Forsyth  (University  of 
Western  Ontario) 

3 

Fremouw  (Physical 
Dynamics,  Inc.  ) 

4 

Rino  (SRI  International) 

Plasma 

Ion  and  electron 

^Langmuir  probe:  electron 

^Rich  (Air  Force 

monitor 

concentration 

concentration  and 

Geophysics  Laboratory) 

Dominant  ion  species 

temperature;  spectral 

Heelis.  Hanson 

Ion  and  electron  temperature 

density  of  concentration 

(University  of 

Electric  field 

irregularities  at  3,  10, 

35,  500,  and  750  m 
wavelengths. 

Retarding  potential  analyzer; 
ion  concentration  and 
in-track  ion  velocity, 

3,  5  km  resolution; 
temperature  and  mass  of 
dominant  ion. 

Ion  drift  meter:  cross-track 
drift  velocity  field  with 

880  m  resolution. 

Texas.  Dallas) 

Vector 

Magnetic  field 

Three -axis  fluxgate:  13  nT 

4Potemra  (A PL) 

magnetometer 

Electric  current  along 
magnetic  field  lines 

quantizations.  400  m 
resolution. 

Electron 

Electron  flux  and  energy 

Energy  range:  20  eV  to 

4 

Hardy  (Air  Force 

spectrometer 

spectra 

20  keV;  three  look 
directions:  up,  down, 
and  45°  from  zenith. 

Three  modes,  spatial 
resolution  300,  600.  and 

1800  m. 

Geophysics  Laboratory) 

Auroral 

Night  and  day  auroral 

2  t 

1150  to  2000  A  imaging 

4 

Huffman  (Air  Force 

Ionospheric 

and  F -layer  emissions 

spectrophotometer. 

Geophysics  Laboratory) 

Mapper 

Wavelength  resolution:  30  A; 
spatial  resolution  in  E  (F) 
layer  is  5  X  20  km 
(3  X  13  mi)  at  nadir. 

3914  and  6300  A  nadir- 
viewing  photometers. 

Meng  (APL) 

'Not  fully  operable.  3H1LAT  Project  Scientist.  L.  A.  Wittwer  (Defense  Nuclear 

_  .  Agency),  Program  Manager 

^Imager  failed  on  23  July  1983.  ^Experiment  principal  investigator. 
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The  major  objective  of  the  HiLat  satellite  program  is  to  provide  a  quantitative 
specification  of  high  latitude  scintillation  strength  and,  in  particular,  the  temporal 
and  spatial  spectra  of  intensity  and  phase  fluctuations  and  the  shape  of  the  irregu¬ 
larity  structures.  These  parameters  provided  by  the  radio  beacon  experiments  are 
supported  by  the  simultaneous  in-situ  data  that  define  the  background  ionospheric 
processes. 

In  this  report  we  shall  concentrate  on  the  HiLat  satellite  observations  per¬ 
formed  by  the  Air  Force  Geophysics  Laboratory  at  Tromso,  Norway  during  the 
1984-1985  period.  The  station  is  located  in  the  central  part  of  the  auroral  oval 
during  the  nighttime  hours  under  magnetically  quiet  conditions  as  defined  by  the 
planetary  magnetic  index  Kp  <  3.  5.  We  shall  study  the  strength  and  the  structure 
specifications  of  complex  signal  scintillations  in  the  spatial  and  temporal  frames 
of  invariant  latitude  and  magnetic  local  time.  These  distributions  will  then  be 
examined  in  the  context  of  in-situ  plasma  structures  observed  by  the  HiLat  satellite 
in  this  environment  and  the  theoretical  predictions  of  structures  that  characteristic¬ 
ally  develop  in  this  region  through  gradient  drift  and  Kelvin-Helmholtz  instability 
processes. 


2.  DATA 

We  shall  illustrate  the  space-time  variations  of  the  statistical  parameters  that 
define  the  complex  scintillation  magnitudes  and  their  structure  over  the  two-year 
period  December  1983  -  October  1985.  Scintillation  measurements  have  been  made 
on  a  sliding  30-sec  window  with  values  computed  every  15  sec  at  the  midpoint  of 
the  window  over  the  length  of  each  pass  (approximately  10  min)  and  have  been 
merged  into  seasonal  data  bases.  To  avoid  the  effects  of  multipath  propagation, 
the  data  acquired  above  a  satellite  elevation  angle  of  20°  were  used  in  the  study. 
This  provided  a  maximum  latitude  coverage  of  about  ±7°  at  an  ionospheric  heigh* 
of  350  km  around  the  station.  Due  to  the  precession  of  the  Hilat  orbit  to  earlier 
times  each  day,  a  full  coverage  in  local  time  is  also  obtained  over  a  season  (that 
is,  three  months)  when  the  ascending  and  descending  node  passes  are  combined. 

The  parameters  discussed  in  this  report  are  the  137-MHz  phase  and  intensity  scin¬ 
tillation  observations  and  their  associated  spectral  parameters,  the  total  electron 
content,  and  the  in-situ  ion  density  measured  at  the  height  of  the  satellite.  These 
statistical  parameters  used  to  quantify  scintillations  are  analogous  to  those  defined 
in  earlier  Wideband  (Basu  et  al ) ^  and  Fleetsat  (Basu  et  al)^  reports.  A  short 
definition  of  each  follows. 
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2.1  RMS  Phase  Deviation 


The  standard  deviation  of  the  linear  detrended  phase  component  over  a  30-sec 
data  interval  is  calculated  as 

a<f)  =  (<*2>  -  <*»2>1/2 
and  is  the  phase  scintillation  parameter. 

22,  Phase  Spectral  Slope  and  Strength 

The  phase  power  spectral  density  ($^(f))  assumes  a  power  law  form  as 

'prf> 

$Af)  =  T.  f  9 

0  0 

where  f  is  the  fluctuation  frequency,  p^  is  the  spectral  slope  and  T^  is  the  spectral 
strength  in  mks  units.  The  slope  is  computed  from  a  linear  least-squares  fit  to 
the  plot  of  the  logarithm  of  4>^(f)  vs  the  logarithm  of  frequency  over  the  range 
0.  2  to  10  Hz.  The  spectrum  is  generated  by  a  FFT  algorithm  with  a  cosine  window 
on  the  30-sec  data  sample. 

23  Inten&'ty  Scintillation  Index, 

The  index  was  first  introduced  by  Briggs  and  Parkin11  and  is  the  normalized 
second  central  moment  of  signal  intensity: 

S4  =  [(<I2>  -  <I>2)/<I>2]1/2. 

It  ranges  from  0  to  ~  1  for  quiet  to  strong  scatter  conditions. 

2.4  Decorrelation  Time 

Decorrelation  time  is  computed  as  the  time  delay  in  seconds  for  the  intensity 
signal  to  reach  a  decorrelation  level  of  0.  5. 


11.  Briggs,  B.  H. ,  and  Parkin,  I.  A.  (1963)  On  the  variation  of  radio  star  and 
satellite  scintillation  with  zenith  angle,  J,  Atmos.  Terr.  Phys.  , 
25:339. 
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The  phase  and  intensity  scintillation  parameters  pertain  to  the  350-km  sub- 

ionospheric  intersection  point  of  the  satellite  as  observed  from  Tromso.  The  total 

electron  content  (TEC)  is  measured  up  to  830  km,  the  height  of  the  satellite,  and 

1/2 

the  ion  density  is  that  measured  at  the  height  of  the  satellite,  830  km.  A  (sec  i) 
correction  for  ionospheric  zenith  angle  (i)  was  incorporated  into  the  a  and 
parameters  while  a  (sec  i)  correction  was  used  for  TEC.  Median  values  of  the  above 
mentioned  parameters  have  been  obtained  in  2.  5°  invariant  latitude  -  1  hr  MLT 
bins  and  will  be  presented  in  polar  plot  form  in  the  following  sections.  A  maximum 
of  150  data  points  are  available  in  some  invariant  latitude-MLT  bins  but  the 
majority  of  the  bins  have  50  data  points.  Bins  with  fewer  than  10  data  samples 
have  been  omitted  to  provide  some  statistical  validity  to  the  results. 

The  invariant  latitude  of  the  overhead  position  of  Tromso  is  66°  A  which  corres¬ 
ponds  to  the  central  part  of  the  nightside  auroral  oval.  Figure  1  shows  the  station 

1 2 

location  at  midnight  with  respect  to  the  Q  =  3  auroral  oval  (Feldstein  and  Starkov), 

13 

which  is  reproduced  from  Whalen.  This  allows  comparison  of  the  day-  and 
nightside  auroral  oval  locations  under  quiet  magnetic  conditions  with  the  band  of 
latitudes  for  which  data  is  available  (usually  60°  -72.  5°  A).  The  invariant  latitudes 
agree  almost  exactly  with  the  corrected  geomagnetic  latitude  (used  for  Figure  1)  in 
this  longitude  sector  (J.  A.  Whalen,  private  communication,  1984).  Thus  we  note 
that  the  HiLat  data  from  Tromso  generally  covers  the  nightside  oval  while  the  day- 
side  oval  is  immediately  poleward  of  the  coverage  from  the  station. 

The  seasonal  data  bases  have  been  sorted  into  two  groups:  those  observations 
made  under  magnetically  quiet  conditions  (Kp  <  3.  5)  and  those  under  magnetically 
disturbed  conditions  (Kp  >  3.5).  The  data  sub-set  under  magnetically  disturbed 
conditions  (Kp  >  3.  5)  is  very  limited.  This  constrains  us  to  select  only  one  season 
as  the  "worst  case"  and  illustrate  the  median  values  for  all  the  parameters  during 
this  season  under  conditions  of  Kp  >  3.  5. 


12.  Feldstein,  Y.  I.  ,  and  Starkov,  G.  V.  (  1967)  Dynamics  of  auroral  belt  and  polar 

geomagnetic  disturbances.  Planet.  Space  Sci.  ,  15:209. 

13.  Whalen,  J.A.  (197  0)  Auroral  Oval  Plotter  and  Nomograph  for  Determining 

Corrected  Geomagnetic  Local  Time,  Latitude  and  Longitude  for  High 
Latitudes  in  Northern  Hemisphere.  AFCRL-70-0422,  AD  713170, 

Air  F"orce  Geophysics  Laboratory,  Hanscom  AFB,  MA. 
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AURORAL  OVAL 
(Feldstein  ond  Storkov,  1967) 


Figure  1.  Polar  Plot  of  Invariant  Latitude  and  Magnetic 
Local  Time  Illustrating  the  Latitude  Interval  of  Available 
Data  With  Respect  to  the  Q=  3  Auroral  Oval 


One  remaining,  but  very  important,  factor  to  be  considered  is  the  solar  flux 
variation  over  the  two-year  period.  Figure  2  illustrates  the  drastic  variation  of 
sunspot  number  for  Dec  1983  -  Oct  1985.  As  a  reference,  at  the  peak  of  the  last 
sunspot  cycle  in  the  fall  of  197  9,  sunspot  numbers  approached  values  of  about 
180  while  in  the  fall  of  1985,  sunspot  numbers  decreased  to  values  as  low  as  20. 


Figure  2.  Variation  of  Monthly  Mean  Sunspot  Number  Between 
December  1983  and  December  1985 


3.  RESULTS 

3.1  RMS  Phase  Deviation 

Figure  3  illustrates  the  rms  phase  deviation,  o^,  of  137  MHz  scintillation  for 
winter  1983  (start  of  observations  in  Dec  1983  to  Jan  1984),  spring  1984  (Feb-Apr 
1984),  summer  1984  (May-Jul  1984)  and  fall  1984  (Aug-Oct  1984)  for  Kp  <  3.5.  The 
rms  phase  deviation  is  computed  over  a  30-sec  data  interval  and  the  variation  of  the 
parameter  with  ionospheric  zenith  angle  has  been  removed  prior  to  statistical 
processing.  Since  the  projected  scan  velocity  of  the  satellite  is  approximately 
3  km/sec  at  350  km,  the  data  interval  covers  irregularity  scale  sizes  as  large  as 
90  km.  The  data  are  binned  in  2.  5°  latitude  and  1  hr  MLT  intervals  and  the  median 
values  are  indicated  for  each  bin.  Each  bin  contains  at  least  ten  data  points.  The 
tic  marks  along  the  noon-midnight  and  the  dawn-dusk  meridians  indicate  10°  intervals 
between  50°  -  90°  invariant  latitudes.  The  sharp  increase  in  rms  phase  deviation 
exceeding  5  radians  over  a  narrow  latitude  swath  (65°  -67.  5°  A)  between  the 
pre -midnight  and  dawn  periods  is  the  most  conspicuous  feature  of  the  plots.  This 
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Figure  3.  Seasonal  Variation  of  Median  rms  Phase  Deviation,  ca  (137  MHz),  With 
Invariant  Latitude  and  Magnetic  Local  Time  Under  Quiet  Magnetic  Conditions 
During  1984 


region  of  enhancement  corresponds  to  the  region  where  the  alignment  of  the  ray 

14 

path  with  the  magnetic  L-shell  occurs  (Fremouw  et  al).  This  location  also  coin- 

15 

cides  with  the  position  of  the  diffuse  aurora  (Hardy  et  al),  a  point  we  will  discuss 

16 

later.  It  is  interesting  to  note  that  in  a  companion  report,  Weimer  has  shown  that 


14.  Fremouw,  E.  J. ,  Rino,  C.L.,  Livingston,  R.  C.  ,  and  Cousins,  M,  D.  (1977) 

A  persistent  subauroral  scintillation  enhancement  observed  in  Alaska, 
Geophys.  Res.  Lett. ,  4:539. 

15.  Hardy,  D.A.,  Gussenhoven,  M.S.,  and  Holeman,  E.  (1985)  A  statistical 

model  of  auroral  electron  precipitation,  J.  Geophys.  Res. ,  90:4229. 

16.  Weimer,  D.  R.  (1987)  Large-scale  Plasrra  Density  Fluctuations  Measured  With 

the  HiLat  Satellite  at  835  km  Altitude,  ATGL-TR-37-01iO  ADA  183043, 

Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA. 
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medium  scale  blobs  (~60  km)  withAN/N  of  10  percent  maximize  at  exactly  the 
same  2.5°  latitude  bin  on  the  nightside.  It  is  possible  that  km-scale  irregularities 
are  associated  with  the  edges  of  these  blobs.  Thus  a  combination  of  high  irregu¬ 
larity  amplitude  and  favorable  geometry  contributes  to  the  pronounced  increase  of 
phase  scintillations.  The  magnitude  of  phase  scintillations  during  winter  is  found 
to  be  larger  than  summer.  This  appears  somewhat  intriguing  because  the  TEC  (to 
be  shown  later)  was  higher  in  summer.  This  leads  us  to  conclude  that  the  irregu¬ 
larities  are  probably  less  preponderant  during  summer  in  the  presence  of  enhanced 
ionization  of  the  underlying  E- region.  Larger  phase  deviations,  o^,  and  hence 
increased  plasma  density  deviation  observed  during  the  vernal  equinox  of  1984  can 
be  related  to  increased  solar  activity  during  this  period  as  will  be  shown  in  a  sub¬ 
sequent  diagram. 

Figure  4  presents  the  rms  phase  deviation,  o^,  obtained  during  1984  and  1985 
in  a  different  format.  This  figure  indicates  the  variation  of  median  rms  phase 
deviation,  o^,  with  invariant  latitude  in  the  midnight  and  noon  time  periods  only. 
Basically,  this  is  obtained  from  a  midnight-noon  cut  through  the  dial  plots  of 
Figure  3.  When  the  number  of  data  points  in  a  particular  bin  over  the  noon-midnight 
sector  falls  below  ten,  the  values  are  extrapolated  from  an  adjacent  time  sector. 
Figure  4  illustrates  an  overall  decrease  of  scintillations  in  1985  due  to  decreased 
solar  activity. 


3.2  Phase  Spectral  Strength. 

Figure  5  shows  the  behavior  of  the  spectral  strength,  T^,  of  phase  scintillations 
during  1984.  The  phase  spectral  strength  is  obtained  from  the  phase  spectra  and 
corresponds  to  the  power  spectral  density  at  1  Hz.  A  midnight-noon  cut  of  the 
spectral  strength  dial  plot  (Figure  5)  is  shown  in  Figure  6,  with  the  addition  of 
1985  data.  The  levels  of  T^  in  the  diagram  may  be  properly  interpreted  by  con¬ 
sidering  that  the  noise  level  of  the  phase  spectra  corresponds  to  about  -50  dB. 

Figure  4  indicates  that  T^  becomes  maximum  in  the  nighttime  sector  over  the  lati¬ 
tude  swath  of  65°  -67.5®  A  where  the  rms  phase  deviation,  a also  becomes  maxi¬ 
mum  as  discussed  in  the  previous  section.  However,  the  enhanced  T^  is  consider¬ 
ably  more  extended  in  latitude  than  a It  seems  that  T^  being  evaluated  at  1  Hz 
corresponds  to  an  irregularity  scale  size  of  about  3  km  whose  sheet-like  anisotropy 
is  less  pronounced  than  the  larger  irregularities  with  tens  of  km  scale  sizes  that 
control  the  magnitude  of 
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Figure  4.  Seasonal  Variation  of  Median  rms 
Phase  Deviation,  (137  MHz),  With 
Invariant  Latitude  ^Along  the  Midnight-Noon 
Meridian  Under  Quiet  Magnetic  Conditions 
in  1984  and  1985 
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Figure  5.  Seasonal  Variation  of  Median  137  MHz  Phase  Spectral  Strength,  T. , 
With  Invariant  Latitude  and  Magnetic  Local  Time  Under  Quiet  Magnetic  ^ 
Conditions  During  1984 
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PHASE  SPECTRAL  STRENGTH 


137MHz  MEDIAN  PHASE  SPECTRAL  STRENGTH  ,  T^  Kp<3.5 

-  1984 


Figure  6.  Seasonal  Variation  of  Median  137  MHz  Phase 
Spectral  Strength,  T  ,  With  Invariant  Latitude  Along 
the  Midnight-Noon  Meridian  Under  Quiet  Magnetic 
Conditions  in  1984  and  1985 
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3.3  Phase  Spectral  Slope,  p0 


Figure  7  shows  the  behavior  of  the  power  law  spectral  index,  p^,  for  phase 
scintillations  at  137  MHz  during  1984.  A  linear  least  square  fit  to  the  phase 
spectrum  indicating  the  variation  of  the  logarithm  of  power  spectral  density  (psd) 
with  the  logarithm  of  frequency  is  obtained  over  the  frequency  interval  of  0.  2  Hz 
to  10  Hz  as  mentioned  earlier.  The  best-fit  line  provides  p ,  as  it  defines  the 

-p .  * 

dependence  of  psd  on  frequency  f  as  psd  cc  f  Considering  the  scan  velocity  of 
the  ray  path  through  the  F- region,  the  fit  range  0.  2  Hz  to  10  Hz  corresponds  to 
the  scale  size  regime  of  about  15  km  to  300  m  for  isotropic  irregularities.  This 
regime  encompasses  the  dominant  structures  that  cause  phase  and  intensity  scin¬ 
tillations  at  VHF  and  UHF  over  the  observing  data  interval.  Figure  8  presents 
the  midnight-noon  cut  of  the  dial  plot  of  power-law  spectral  index,  p^,  with  the 
addition  of  the  1985  data. 

From  Figures  7  and  8,  it  may  be  noted  that  the  phase  spectral  index  is  least 

affected  by  the  propagation  geometry  and  does  not  show  significant  variations  with 

season  and  solar  activity  in  contrast  to  the  behavior  of  o  , .  There  exists,  however, 

<P 

a  tendency  for  p^  to  increase  around  65°  A,  during  the  nighttime.  Since  this  loca¬ 
tion  corresponds  to  the  average  location  of  the  diffuse  auroral  oval  as  mentioned 
earlier,  the  associated  E- region  conductivity  may  account  for  the  damping  of  the 

short  scale  irregularities  and  causing  the  phase  spectrum  to  steepen  (Vickrey  and 
17 

Kelley).  This  steepening  may  also  arise  from  increased  psd  at  larger  scales  in 
the  geometrical  enhancement  region  due  to  better  L-shell  alignment  of  large  scale 
irregularities  as  noted  in  connection  with  the  increased  o^.  The  enhancements  of 
p^  are,  however,  smooth  and  distributed  in  contrast  to  the  sharp  geometrical 
enhancements  of  a Probably,  the  spectral  steepening  arises  from  an  interplay 
of  both  geophysical  and  geometrical  effects. 


17.  Vickrey,  J.F.,  and  Kelley,  M.C.  (1982)  The  effects  of  a  conducting  E-layer 
on  classical  F -region  cross-field  plasma  diffusion,  J.  Geophys.  Res.  , 
87:4461. 
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Figure  7.  Seasonal  Variation  of  Median  137  MHz  Phase  Spectral  Slope, 
pu  (0.  2  -  10  Hz),  With  Invariant  Latitude  and  Magnetic  Local  Time  Under 
Qhiet  Magnetic  Conditions  During  1984 
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Figure  8.  Seasonal  Variation  of  Median  137  MHz 
Phase  Spectral  Slope,  pu  (0.  2-  10  Hz),  With 
Invariant  Latitude  Along  the  Midnight-Noon 
Meridian  Under  Quiet  Magnetic  Conditions  in  1984 
and  1985 
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3.4  Intensity  Scintillation  Index,  S4 

The  intensity  scintillation  magnitudes  have  been  expressed  in  terms  of  the 

standard  S.  index,  defined  as  the  normalized  variance  of  signal  intensity  (Briggs 
^11  1/2 

and  Parkin).  As  mentioned  earlier,  a  (sec  i)  variation  of  with  ionospheric 
zenith  angle  has  been  removed  prior  to  statistical  processing.  Figure  9  shows 
the  variation  of  median  at  137  MHz  with  invariant  latitude  and  MLT.  The  night¬ 
time  enhancement  of  intensity  scintillation  index,  S^,  in  the  region  of  alignment  of 
the  ray  path  with  the  magnetic  L-shell  (65°  -  67.  5°A)  is  observed  to  be  much  less 
pronounced  when  it  is  compared  to  the  rms  phase  deviation  o^,  enhancements 
shown  in  Figure  3.  This  indicates  that  the  L-shell  alignment  of  km-scale  irregu¬ 
larities  causing  intensity  scintillations  is  considerably  less  than  the  irregularities 
in  the  tens  of  kilometers  scale  that  control  the  o.  values. 

Figure  10  shows  the  midnight-noon  cut  of  the  dial  plots  of  Figure  9  with  the 
addition  of  1985  data.  The  decrease  in  intensity  scintillation  with  decrease  in 
solar  activity  is  in  keeping  with  that  of  the  rms  phase  deviation  shown  in  Figure  4. 

Putting  these  intensity  scintillation  index,  S.,  results  in  perspective  with  other 

3  4 

data  sets,  Basu  and  Aarons  had  hypothesized  that  the  seasonal  variation  of  in 
the  Scandinavian  sector  would  not  be  as  dramatic  as  seen  at  Narssarssuaq, 

Greenland  (64°  A).  The  Narssarssuaq  data  base,  compiled  from  137  MHz  observa¬ 
tions  of  ATS-3  during  1968-1974  (a  high  solar  flux  period),  showed  a  strong  seasonal 
pattern.  The  ATS-3  observations  indicated  that  during  the  Feb  -  Apr  period  a  night¬ 
time  maximum  value  of  mean  S.  =  0.  5  was  obtained  which  decreased  to  S.  =  0.  25  in 

4  18  ^ 
winter.  More  recently,  Kersley  et  al  have  shown  observations  from  Kiruna  at 

150  MHz  of  intensity  scintillation,  S^,  during  Sept  1984  -  Dec  1985.  In  this  data 

base,  the  nighttime  maximum  over  the  range  60  -  72.5°  CGMLAT  occurs  at  ~65° 

with  median  S4  -  0.2.  The  winter  data  shows  a  slightly  lower  median  S^. 


18.  Kersley,  L. ,  Pryse,  S.E.,  and  Wheadon,  N.  S.  (1986)  Radio-wave 

Scintillations  and  Ionospheric  Irregularities  at  High  Latitudes,  Rpt.  Depart¬ 
ment  of  Physics,  University  College  of  Wales,  Aberystwyth,  U.  K. 
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Figure  9.  Seasonal  Variation  of  Median  Intensity  Scintillation  Index,  S4,  at 
137  MHz  With  Invariant  Latitude  and  Magnetic  Local  Time  Under  Quiet 
Magnetic  Conditions  During  1984 
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Figure  10.  Seasonal  Variation  of  Median  Intensity 
Scintillation  Index,  S4,  at  137  MHz  With  Invariant 
Latitude  Along  the  Midnight- Noon  Meridian  Under 
Quiet  Magnetic  Conditions  in  1984  and  1985 
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3.5  Decorrelation  Time,  r 


Decorrelation  time,  r,  is  computed  as  the  time  delay  in  seconds  for  the  inten¬ 
sity  signal  to  reach  a  decorrelation  level  of  0.5.  For  orbiting  satellites,  such  as 
HiLat,  the  decorrelation  time  is  dictated  by  the  effective  velocity  of  the  satellite 
at  F -region  heights  as  well  as  the  electron  density  structure.  The  effective  velocity 
is  determined  by  the  irregularity  anisotropy  and  the  scan  velocity  of  the  ray  path 
through  the  ionosphere  (Rino).  The  midnight  cut  of  decorrelation  time  for  the 
available  range  of  invariant  latitudes  is  shown  in  Figure  11.  When  the  intensity 
scintillation  index  is  low,  decorrelation  times  cannot  be  measured  accurately. 
Hence,  daytime  decorrelation  times  are  not  shown.  The  decrease  of  decorrelation 
time  to  values  as  low  as  0.  25  is  obtained  at  65  -  67.  5°  A  where  the  intensity  scintilla¬ 
tion  index,  S^,  becomes  maximum.  This  position  corresponds  to  the  alignment  of 
the  ray  path  the  magnetic  L-shell.  The  effective  velocity  becomes  maximum  at  the 
above  position  which  probably  causes  a  decrease  of  the  decorrelation  time. 

3.6  TEC  (Total  Electron  Content) 

The  comb  of  three  UHF  transmissions  from  HiLat  is  used  to  derive  the  total 
electron  content  of  the  ionosphere  up  to  the  satellite  altitude  of  830  km  by  the  differ- 
ential  Doppler  technique  (Fremouw  et  al).  In  Figure  12  is  shown  the  variation  of 

median  TEC  with  invariant  latitude  and  MLT  while  Figure  13  shows  the  latitude 
variation  of  median  TEC  values  during  the  noon  and  midnight  cut  of  the  dial  plot 
shown  in  Figure  12  for  different  seasons  in  1984  and  1985.  The  next  diagram.  Fig¬ 
ure  14,  shows  the  observed  variation  of  noontime  TEC  in  relation  to  the  sunspot 
number.  The  effects  of  solar  activity  and  season  on  TEC  appear  to  be  coupled.  At 
high  latitudes,  owing  to  the  near  vertical  orientation  of  the  earth's  magnetic  field, 
the  ionospheric  irregularities  of  electron  density  at  F-region  heights  usually  extend 
in  altitude.  The  integrated  effect  of  the  irregularities  on  radio  wave  propagation, 
such  as  scintillation,  is  therefore  weighted  by  the  total  electron  content  of  the  ionos¬ 
phere.  Thus  measurements  of  irregularity  amplitude,  AN/N  (AN  being  the  rms 
electron  density  fluctuation  and  N  the  background  density),  and  total  electron  content 
may  form  the  basis  for  a  modeling  of  scintillation  magnitudes. 

Median  foF2  data  from  Sodankyla  (63°  CGM)  (S.  Kirkwood,  private  communica¬ 
tion,  1986)  shows  a  similar  nighttime  pattern  with  Dec  1983  -  Jan  1984  and 


19.  Rino,  C.L.  (1979)  A  power  law  phase  screen  model  for  ionospheric  scintillation. 

1,  Weak  scatter.  Radio  Sci. ,  14: 1 135. 

20.  Fremouw,  E.  J. ,  Leadabrand,  R.  L. ,  Livingston,  R.C.,  Cousins,  M.  D.  , 

Rino,  C.L.,  Fair,  B.C.,  and  Long,  R.A.  (1978)  Early  results  from  the 
DNA  Wideband  satellite  experiment  -  complex- signal  scintillation. 

Radio  Sci. ,  13:167. 


Nov  1984  -  Jan  1985  showing  minimum  values  (~  2 -3  MHz)  while  May-.lul  1984 
is  the  period  of  maximum  values  (~  5  MHz).  Daytime  foF2  shows  much  less  of  a 
seasonal  pattern  overall  with  maximum  in  Feb  -  Apr  1984  at  8  MHz. 


137MHz  MEDIAN  DECORRELATION  TIME  Kp<3.5 


-  1984 

OOMLT  - 1985 


Figure  11.  Seasonal  Variation  of  Median  Decorrelation  Time,  r  , 
(137  MHz)  With  Invariant  Latitude  Along  the  Midnight-Noon  Meridian 
Under  Quiet  Magnetic  Conditions  in  1984  and  1985 
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Figure  12.  Seasonal  Variation  of  Median  Total  Electron  Content  (TEC)  With 
Invariant  Latitude  and  Magnetic  Local  Time  Under  Quiet  Magnetic  Conditions 
During  1984 


22 


MEDIAN  TEC 


Kp  <3.5 


-  1984 

- 1985 


OOMLT 


12  MLT 


INVARIANT  LATITUDE 


Figure  13.  Seasonal  Variation  of  Median  Total 
Electron  Content  (TEC)  With  Invariant  Latitude 
Along  the  Midnight-Noon  Meridian  Under  Quiet 
Magnetic  Conditions  in  1984  and  1985 
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Figure  14.  Illustrates  the  Correspondence  Between  the  Variations 
of  Total  Electron  Content  and  Sunspot  Number 
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3.7  Ion  Density 


Among  the  various  in-situ  parameters  probed  by  the  satellite -born  sensors, 
the  in-situ  ion-density  measurement  is  an  important  parameter.  In  the  F-region, 
due  to  charge  neutrality,  the  electron  and  ion  densities  are  equal.  Figure  15 
illustrates  the  invariant  latitude  -  MLT  variation  of  ion  density  during  1984  while 
Figure  16  shows  the  variation  of  this  parameter  with  latitude  in  the  midnight-noon 
time  frame  as  a  function  of  season  in  both  1984  and  1985  under  quiet  magnetic 
conditions.  The  latitudinal  coverage  in  both  Figures  15  and  16  is  wider  than  seen 
in  the  previous  figures  since  the  ion  density  measurement  is  an  in-situ  parameter. 
The  daytime  ion  density  at  830  km  is  observed  to  follow  closely  the  pattern  of  TEC 
variations  shown  earlier  in  Figures  12  and  13.  This  is  seen  in  Table  2. 


Table  2.  TEC  and  Ion  Density  Measurements 
67.5°  A  12  MLT 


TEC 

<x  10"16  m‘2) 

Ion  Density 

,  in-10  -3, 

(x  10  m  ) 

Nov  -  Jan  1984 

5 

0.6 

Feb  -Apr  1984 

19 

3.3 

May  -Jul  1984 

12 

2.3 

Aug  -  Oct  1984 

3 

1.  o 

Nov  -  Jan  1985 

1 

0.3 

Feb  -  Apr  1985 

8 

1.  0 

May  -  Jul  1985 

9 

1.  3 

Aug  -Oct  1985 

3 

0.  6 

The  ion  density  variations  at  830  km  are  observed  to  vary  with  both  solar  activity 
and  season. 
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Figure  15.  Seasonal  Variation  of  Median  Ion  Density  With  Invariant  Latitude 
and  Magnetic  Local  Time  Under  Quiet  Magnetic  Conditions  During  1984 
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Table  3  (reproduced  from  Rich  and  Smiddy)  illustrates  a  comparison  of  ion 
density  observed  under  solar  maximum  conditions  using  the  DMSP  satellite  in  1979 
with  the  corresponding  ion  density  observed  under  solar  minimum  conditions  in 
1984  with  the  HiLat  satellite. 


Table  3.  Comparison  of  Median  Density  at  60°  Magnetic  Latitude 
Observed  at  830  km  Altitude  by  DMSP  in  1979  and  HiLat  in  1984 
(Reproduced  from  Rich  and  Smiddy)^1 


In-situ  Density 

in-4  ~3 

x  10  cm 

Summer 

1979 

DMSP 

Winter 

1984  1979 

HiLat  DMSP 

1984 

HiLat 

Morning 

07  -  11  hrs  MLT 

10 

1.8  2.5 

0.  18 

Evening 

19  -  23  hrs  MLT 

10 

2.  5  4  to  7 

0.36 

33  Effect  of  Magnetic  Activity 

Magnetic  activity,  in  addition  to  sunspot  number,  exerts  overwhelming  control 

3  9 

of  scintillations  in  the  auroral  region  (Basu  and  Aarons  and  Basu  et  al  ).  The 
"worst  case"  conditions  can  be  illustrated  by  the  Feb  -  Apr  1984  season  when  sun¬ 
spot  numbers  higher  than  those  seen  in  the  remainder  of  1984  and  1985  are  coupled 
with  magnetically  disturbed  conditions. 

Our  HiLat  observations  indicate  that,  under  magnetically  active  conditions, 
scintillation  magnitudes  are  enhanced  during  the  nighttime  and  the  region  of  en¬ 
hanced  activity  extends  in  both  the  poleward  and  the  equatorward  directions.  Ir- 
3spective  of  the  level  of  solar  activity,  scintillations  are  enhanced  during  mag- 
tic  dL  irbances.  In  Figure  17  is  shown  the  enhancement  of  the  following 
..nr  v.  ■  rs  under  magnetically  disturbed  conditions  during  Feb- Apr  1984:  rms 
phase  deviation,  phase  spectral  strength,  T  ;  phase  spectral  slope,  p^; 
intensity  scintillation  index,  S^;  and  decorrelation  time. 

"  1 .  Kich,  F.J. ,  and  Smid  ,  M.  (1986)  Plasma  Densities  and  Irregularities  at 

830-km  Altitude  Bat  d  on  Observations  During  1979,  AFGL  Technical  Rpt. 
^FGL-TH-66-0l2i,  ADA  172115,  Air  Force  Geophysics  Laboratory, 

!  ’scorn  AFB,  MA. 
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137MHz  FEB -APR  1984  MEDIAN  VALUES  Kp>3.5 


Figure  17.  Variations  of  Scintillation  Parameters:  rms  Phase  Deviation,  a 
Phase  Spectral  Strength,  Ta;  Phase  Spectral  Slope,  p^  Intensity  Scintillation 
Index,  S4;  and  Decorrelation  Time,  r.  Under  "Worst  Case"  Conditions. 

These  observations,  made  during  Feb  -  Apr  1984  under  disturbed  magnetic 
conditions,  are  shown  as  a  function  of  invariant  latitude  along  the  midnight-noon 
meridian 
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Under  the  "worst  case"  conditions,  both  the  rms  phase  deviations,  and 

the  intensity  scintillation  index,  S^,  increase  by  50  percent  over  the  scintillation 

level  observed  under  quiet  magnetic  conditions  for  Feb -Apr  1984.  The  phase 

spectral  slope,  p^,  becomes  steeper  and  the  phase  spectral  strength,  T^,  greater 

under  disturbed  magnetic  conditions.  No  change  is  seen  in  decorrelation  time,  t, 

possibly  due  to  the  fact  that  the  effective  scan  velocity  controls  t  for  all  values  of 

S.  <  0.  5. 

4 

The  large  increase  in  cr^,  T^  and  in  the  70°  -  72.  5°  A  latitude  region  in  the 
daytime  during  magnetic  disturbances  is  in  all  probability  caused  by  the  equator- 
ward  motion  of  the  dayside  cusp  compared  to  its  location  shown  in  Figure  1. 

4.  DISCUSSION 

We  have  shown  that  the  most  notable  feature  in  the  observed  scintillation 
morphology  at  this  auroral  station  even  under  magnetically  quiet  to  moderate  con¬ 
ditions  is  the  enhancement  of  scintillations  over  a  narrow  latitude  interval  centered 
at  65°  A  in  the  nighttime  sector.  This  narrow  band  increase  of  scintillations  coin¬ 
cides  with  the  expected  region  of  geometrical  enhancement  due  to  the  alignment  of 
the  ray  path  with  the  magnetic  L- shell  oriented  irregularities.  It  is  interesting  to 
note  that  the  boundary  blobs  that  signify  plasma  density  enhancements,  tens  of 

kilometers  wide  in  the  magnetic  E-W  direction,  also  occur  in  this  region  (Rino 
22 

et  al).  The  spatial  configuration  of  the  blobs  is  controlled  by  the  high  latitude 

convection  pattern  that  is  mostly  E-W  in  the  auroral  region  (Heelis  et  al).  *  0 
24 

Robinson  et  al  have  shown,  with  the  help  of  simulation  studies,  that  even  a 
primarily  circular  patch  of  ionization  being  convected  in  from  the  polar  cap  would 
assume  a  narrow  in  latitude  and  elongated  in  longitude  shape.  These  blobs  can 
also  develop  small  scale  irregularities  in  their  trailing  edges  through  EX  B 
instability  mechanisms.  In  general,  the  plasma  blobs  transported  from  distant 
regions  get  continually  structured  in  the  convection  field  and  control  a  major 
source  of  nighttime  auroral  scintillations. 


22.  Rino,  C.L.,  Livingston,  R.  C. ,  Tsunoda,  R.  T. ,  Robinson,  R.  M. ,  Vickrey,  J.F., 

Senior,  C. ,  Cousins,  M.D.,  Owen,  J.,  and  Klobuchar,  J.A.  (1983)  Recent 
studies  of  the  structure  and  morphology  of  auroral  zone  F-region  irregu¬ 
larities,  Radio  Sci. ,  18:1167. 

23.  Heelis,  R.A.,  Lowell,  J.K.,  and  Spiro,  R.W.  ( 1982)  A  model  of  the  high- 

latitude  ionospheric  convection  pattern,  J.  Geophys.  Res.  ,  87:6339. 

24.  Robinson,  R.  M. ,  Tsunoda,  R.  T. ,  Vickrey,  J.F.,  and  Guerin,  L.  (1985) 

Source  of  F-region  ionization  enhancements  in  the  nighttime  auroral  zone, 

J.  Geophys.  Res. ,  90:7533. 
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The  other  notable  finding  related  to  the  scintillation  structure  is  the  steep 

phase  spectral  slopes  observed  in  the  region  of  geometrical  enhancement.  This 

region  coincides  with  the  central  part  of  the  diffuse  aurora  under  magnetically  quiet 

15 

conditions  as  has  been  recently  shown  by  Hardy  et  al  from  DM3P  satellite  observa¬ 
tions.  Unlike  the  limited  HiLat  satellite  data  base  of  particle  precipitation,  the 
DMSP  satellites  have  provided  a  very  large  amount  of  data  which  have  been 
organized  by  Hardy  et  al  to  show  the  characteristic  variations  of  particle  precipita¬ 
tion  with  magnetic  activity.  The  average  energy  of  the  particle  precipitation  in  the 
diffuse  auroral  region  is  on  the  order  of  a  few  keV  that  is  sufficient  to  produce  a 
conducting  E -region.  This  could  reduce  the  lifetime  of  small  scale  irregularities 
(~  100  m)  and  contribute  to  the  observed  steepening  of  the  phase  spectral  indices. 

Basu  et  al  have  illustrated  such  steepening  of  in-situ  density  spectra  in  regions 
of  energetic  auroral  particle  precipitation. 

The  structures  present  in  the  in-situ  ion  density  probed  by  the  HiLat  satellite 

1 6 

have  recently  been  investigated  (Weimer).  This  study  indicates  that  plasma  den¬ 
sity  blobs  at  830  km  are  concentrated  over  a  fairly  narrow  latitude  interval.  This 
region  when  mapped  down  the  magnetic  field  line  corresponds  approximately  with 
the  region  of  enhanced  scintillations.  These  results  indicate  that,  on  a  statistical 
basis,  the  blobs  in  addition  to  their  E-W  extent  are  also  extended  along  the  mag¬ 
netic  field  lines  to  the  topside  ionosphere.  The  altitude  profiles  of  the  blobs  have 

been  published  earlier  using  incoherent  scatter  radar  measurements  (Vick  rev  et 
26 

al). 

The  extension  of  blobs  to  altitudes  as  high  as  830  km  indicates  that  magneto- 

spheric  coupling  effects  need  to  be  included  in  the  analysis  of  high  latitude  F -layer 

instabilities.  The  inclusion  of  ion-inertia  in  the  EXE  instability  theories  results 

in  a  reduction  of  growth  rates  and  generation  of  irregularities  with  spectral  isotropy 
27 

(Mitchell  et  al).  In  the  non-inertial  domain,  the  irregularity  spectra  become 
anisotropic  in  the  north-south  and  the  east-west  directions.  Thus,  the  irregularity 
spectral  indices  are  expected  to  be  different  in  the  north-south  and  east-west 
directions  in  the  non-inertial  case  and  identical  for  the  inertial  case.  At  Tromso, 
in  addition  to  making  HiLat  observations  that  provide  a  north-south  scan  through 
the  ionosphere,  we  have  performed  scintillation  observations  with  near-stationary 


25.  Basu,  Su. ,  Basu,  S. ,  MacKenzie,  E.  ,  Coley,  W.R.  ,  Hanson,  W.  B.  ,  and 

Lin,  C.S.  (1984)  F-region  electron  density  irregularity  spectra  near 
auroral  acceleration  and  shear  regions,  J.  Geophys.  Rea.  ,  89:5554. 

26.  Vickrey,  J.F.,  Rino,  C.  L. ,  and  Potemra,  T.  A.  (1980)  Chatanika/ TRIAD 

observations  of  unstable  ionization  enhancements  in  the  auroral  F-region, 
Geophys.  Res.  Lett.  ,  7:789. 

27.  Mitchell,  H.  G. ,  Jr.,  Fedder,  J.  A. ,  Keskinen,  M.  J.,  and  Zalesak,  S.  T. 

(1985)  A  simulation  of  high  latitude  F  lay«r  instabilities  in  the  presence  of 
magnetosphere  -  ionosphere  coupling,  Geophys.  Res.  Lett.,  12:283. 
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polar  beacon  satellites  that  provide  an  east-west  scan  controlled  by  ionospheric 
motion.  The  spectral  indices  of  both  sets  of  measurements  are  found  to  be 
approximately  equal  indicating  that  the  inertial  effects  are  important  in  the  genera¬ 
tion  of  kilometer  scale  irregularities  at  high  latitudes  that  cause  VHF  scintillations. 

Another  type  of  auroral  irregularity  not  associated  with  large-scale  organized 
density  gradients  (such  as  to  be  found  on  the  edges  of  blobs)  but  associated  with 
velocity  shears  with  shear  gradient  scale  lengths  ~10  km  has  been  identified  from 

HiLat  observations  (Basu  et  al).  These  irregularities  have  somewhat  shallow 
2  Q 

spectra  (Basu  et  al)  as  they  have  considerable  psd  at  the  shorter  scales  and  can 

cause  intense  VHF  scintillations  at  the  edges  of  auroral  arcs  in  association  with 

1  fi 

upward  field-aligned  currents.  Weimer  has  shown  that  velocity  fluctuations 

statistically  maximize  poleward  of  the  medium  scale  blob  structures.  This  is  in 

agreement  with  the  finding  that  the  discrete  aurora  is  found  to  map  to  the  boundary 

plasma  sheet  and  the  diffuse  aurora  (with  which  the  blobs  are  co-located)  maps  to 

30 

the  central  plasma  sheet  in  the  magnetophere  (Winningham  et  al).  It  is  tempting 
to  associate  some  of  the  scintillations  seen  poleward  of  the  geometrical  enhance¬ 
ment  region  and  having  shallower  phase  spectral  slopes  (cf.  Figures  4  and  8)  with 

3  1 

velocity  shears.  Keskinen  et  al  have  recently  shown  that  these  velocity  shears 
can  drive  the  Kelvin- Helmholtz  instability  in  both  collisional  and  collisionless 
environments. 


28.  Basu,  Su. ,  Basu,  S. ,  Senior,  C. ,  Weimer,  D. ,  Nielsen,  E. ,  and 

Fougere,  P.F.  (1986)  Velocity  shears  and  sub-km  scale  irregularities  in 
the  nighttime  auroral  F -region,  Geophys.  Res.  Lett. ,  13:101. 

29.  Basu,  Su. ,  Basu,  S. ,  MacKenzie,  E.  ,  Fougere,  P.  F.  ,  Coley,  W.  R.  , 

Maynard,  N.C.,  Winningham,  J.D.,  Sugiura,  M.  ,  Hanson,  W.  B. ,  and 
Hoegy,  W.  R.  (1988)  Simultaneous  density  and  electric  field  fluctuation 
spectra  associated  with  velocity  shears  in  the  auroral  oval, 

J.  Geophys.  Res. ,  93:115. 

30.  Winningham,  J.  D. ,  Yasuhara,  F. ,  Akasofu,  S. -I. ,  and  Heikkila,  W.  J.  (1975) 

The  latitudinal  morphology  of  10-eV  to  10-keV  electron  fluxes  during 
magnetically  quiet  and  disturbed  times  in  the  2100-0300  MLT  sector, 

J.  Geophys.  Res.  ,  80:3148. 

31.  Keskinen,  M.  J. ,  Mitchell,  H.  G.,  Fedder,  J.A.,  Satyanarayana,  P.  , 

Zalesak,  S.  T. ,  and  Huba,  J.  D.  (1988)  Nonlinear  evolution  of  the 
Kelvin- Helmholtz  instability  in  the  high  latitude  ionosphere, 
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